Optical Spectral Singularities as Threshold Resonances 
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Spectral singularities are among generic mathematical features of complex scattering potentials. 
Physically they correspond to scattering states that behave like zero-width resonances. For a simple 
optical system, we show that a spectral singularity appears whenever the gain coefficient coincides 
with its threshold value and other parameters of the system are selected properly. We explore a 
concrete realization of spectral singularities for a typical semiconductor gain medium and propose 
a method of constructing a tunable laser that operates at threshold gain. 
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Consider an infinite planar slab gain medium that is 
aligned along the cc-axis and the electromagnetic wave 
given by E{z,t) = E e^'^^^~^^^ ex, where is a constant, 
^ is the propagation constant, and stands for the unit 
vector pointing along the positive x-axis. It is easy to 
show that while traveling through the gain medium the 
wave is amplified by a factor of e^^ , where L is the width 
of the gain medium and g is the gain coefficient. The lat- 
ter is related to the imaginary part n of the the complex 
refractive index of the medium. 



n ^ rj + in, 

and the wavelength A := 27rc/a; according to [l|: 
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Often, one places the gain medium between two mirrors 
to produce a (Fabry-Perot) resonator. This extends the 
length of the path of the wave through the gain medium 
and yields a much larger amplification of the wave for the 
resonance frequencies of the resonator. 

In 0, ^6 have outlined an alternative amplifica- 
tion effect that does not involve mirrors. For the system 
we consider, it requires adjusting L and g so that the 
system supports a spectral singularity. This is a generic 
mathematical feature of complex scattering potentials Q 
that obstructs the completeness of the eigenfunctions of 
the corresponding non-Hermitian Hamiltonian operator. 
Physically, a spectral singularity is the energy of a scat- 
tering state that behaves exactly like a zero-width reso- 
nance 0, In this letter, we first reveal the relation- 
ship between spectral singularities and the well-known 
laser threshold condition [l], and then explore the possi- 
bility of tuning the wavelength of the spectral singularity 
by adjusting the pump intensity. This turns the system 
into a tunable laser that operates at the threshold gain. 

It is easy to show that 

E{z, t) = Ee~''^^ij{z)ex, B{z, t) = -ia;~i£;e~'"*V'(^)ey, 



positive y-axis, t/i is a continuously differentiable solution 
of the time-independent Schrodinger equation: 

-TP"{z) + v{z)^{z)^e^{z), (3) 

k := cj/c ~ 27r/A, and v is the complex barrier potential: 



fc^3 for 



\z\ < k. 



Solving the Schrodinger equation ([3]) we find 



B_e-^^'- for z < -L/2 
Q^^^tnkz \z\<L/2 

B+e-''''- for z > L/2, 



(4) 



(5) 



where k S R"*", A- and B- are free complex coefficients 
and Aq, Bq, and i?+ are complex coefficients related 



to A_ and B^. For example, 

^ — iLk r/„ Lk 
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Ml 1, where 



) 2i(n^ - 1) sin(nLfc) 



-2i(n^ - l)sin(nLfc) 
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is the transfer matrix, and for all zi, Z2 G C, 

/(zi,Z2) := e-^"^^^{l + z,Y - - z,f 



(6) 



Because v is an even function of z, the left and right 
reflection and transmission amplitudes coincide. They 
are respectively OTven by i? = — M2i/Af22 = M12/M22 
and T= I/M22, m- 

The spectral singularities are the fc^ values for which 
M22 = 0, 0, i.e., the real k values satisfying 



/(n,M)=0. 



(7) 



Because / is a complex- valued function, Eq. ([7]) is equiv- 
alent to a pair of coupled real transcendental equations 
for three unknown real variables Re(3), Im(3) and afc. In 
Ref. 31 1 outline a method of decoupling these equa- 
tions. Here we give a more direct solution that reveals 
some previously unknown aspects of the problem. 
First, we use ([6]) to express (O in the form 



is a solution of Maxwell's equations for the above system 
provided that iy stands for the unit vector along the 



^ — 2inLk 



1 - n 



= 0. 



(8) 
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Noting that the boundaries of the gain region have a 
reflectivity of 



7^ := 



1 



n+1 



(9) 



we can write © as e'^^^'^" = 7^. Next, we substitute (H]) 
in this equation and use ([2]) to express its left-hand side 
in terms of g. This gives 



(10) 



Taking the modulus of both sides of this equation yields 



-2gL 



m 



(11) 



This is precisely the laser threshold condition [l]. In 
other words, a spectral singularity and the associated 
zero-width resonance appear at the threshold gain: 



9 = 5th 



— In - — —. 
2L TZ^ 



(12) 



We wish to stress that the threshold condition pT|) 
is only a necessary condition for having a spectral sin- 
gularity. It is by no means sufficient. A necessary and 
sufficient condition is Eq. ^ that we can express as 



kL = — —hiTZ. 
2m 



(13) 



A key observation that reveals the discrete nature of spec- 
tral singularities is that, as a complex- valued function, 
\nTZ has infinitely many values; in view of ([T]) and (|9]), 



ln7^ = ln|7^| + 2i 



tan 



2k 



k2 - 1 



(14) 

where m is an arbitrary integer. Using ([l} and ([9]), we 
also find 



\n\ = 



(77- 1)^ 

(rj + 1)2 + k2 



(15) 



Because t] > 0, this equation implies \TZ\ < 1. Further- 
more, in view of (1121). we have 



9 ^ 9th 



(?7-H)2-fK2 
(r/ - 1)2 -I- k2 



(16) 



Next, we return to Eq. (IT^ . Because kL is real, the 
imaginary part of the right-hand side of this equation 
must vanish. This gives 



f]\n\n\ + 2K 



tan 



2k 



7]^ + — 1 

Furthermore, we can express (|13p as 

ln|7^| 



kL 



2k 



0. (17) 



(18) 



Because kL > and \7l\ < 1, HI]) implies k < 0. Ac- 
cording to ([2|) , this corresponds to the situation that the 
medium has a positive gain coefficient. This is a remark- 
able manifestation of the conservation of energy, because 
whenever we arrange the parameters of the system so 
that a spectral singularity is generated, the system be- 
gins emitting radiation. This can happen only for a gain 
medium, i.e., when g > 0. In [3|, we could only demon- 
strate this graphically. Here we have derived it rigorously. 

Eq. (|17|) determines the location of spectral singulari- 
ties in the 77- k plane. In view of the inequalities: 77 > 0, 
K < 0, \TZ\ < 1, and the fact that tan~^ is an odd func- 
tion taking values in (— f , f ), we can satisfy (ITTl) only 
for 771 > 0. It is also instructive to note that solving 
for In 1 7?, I in (fT7|) . substituting the result in (fT8|) . and us- 
ing k — 2'Kv I we find the following expression for the 
frequency v of the spectral singularity. 



tan 



2k 



k' - 1 



(19) 



where := The first term on the right-hand side 
of (jl9p is the usual resonance frequency of a resonator 
of length L with perfectly reflecting boundaries. It is 
important to note that because k and 77 are frequency- 
dependent quantities, there will be certain frequencies for 
which they satisfy ([T7|. A spectral singularity will arise 
if and only if at least one of these frequencies coincides 
with one of the values fulfilling p^ . It turns out that 
if we fix all the physical parameters of the system this 
can happen only for a single critical frequency, i.e., a 
particular mode number m. 

In order to explore this phenomenon suppose that the 
gain medium is obtained by doping a host medium of 
refraction index tiq and that it is modeled by a two- 
level atomic system with lower and upper level popu- 
lation densities Ni and N^, resonance frequency ujq, and 
damping coefficient 7. Then its permittivity (e :— e^n^) 
is given by 



£0 



UJ 



(20) 



where Eq is the permittivity of the vacuum, oj^ := [Ni — 
Nu)e'^ / (meo), and e and m are electron's charge and 
mass, respectively [ll,|l]. In view of ([J), and (PO]). we 
can express ujp in terms of the gain coefficient go at the 
resonance frequency loq. Introducing Aq := 27rc/ajo and 



Ko 



and using 



s/eo 



Apffo 
47r ' 



(77 + ik)^ 



(21) 



(22) 



and ([2]), we have 



w2 2-iKQyJnl + kI cgol 
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FIG. 1: (Color online) Plots of the curves of spectral sin- 
gularities Cm (full and dashed thick purple curves) and the 
dispersion relation C (thin blue curve) for a gain medium 
with specifications (|24|l . The spectral label m of the Cm 
that are displayed ranges over 1400, 1500, 1600, • • • , 2500 for 
the thick full (purple) curves from bottom to top, and 
3000, 4000, 5000, • ■ • , 10000 for the thick dashed (purple) 
curves. As one increases m these curves approach to the 
(5r;-axis from below. The thin dashed (gray) lines are the 
coordinate a^ces. 



For example for a semiconductor gain medium 1] with: 



no 



3.4, Ao = 1500 nm. 



■ 0.02, go = 40cm~\ (24) 



we have kq ^ -4.7747 x 10^* and = -6.4935 x 10"^ 
Substituting these values in (pn|) and using (P^ . we can 
express Srj := rj — no and k as functions of w and use them 
to plot a parametric curve C representing the dispersion 
relation i^U^i. See Figure [H Next, we return to Eq. p7)) 
and consider it as an implicit function for each value of 
m. Figure [1] also shows the graph of these functions for 
several values of m. These are the curves of spectral sin- 
gularities that we label by Cm- As seen from Figured] for 
sufficiently large values of m (i.e., m > 1374), Cm inter- 
sects C at two points. One of these {pm+) has a positive 
St] coordinate and corresponds to a frequency LUm+ that 
is larger than ujq. The other {pm+) has a negative Sr] and 
corresponds to a frequency Wm- that is smaller than luq. 
We can compute the coordinates (5rim±,Km±) of Pm±^ 
find ujm± and km± '■= ^m±lc^ and use to determine 
L. This gives rise to a set of values Lm± for the length 
of the gain medium that would allow for the existence of 
a spectral singularity. For the gain medium ([M)) , setting 
m = 1374 we find a pair of spectral singularities with 
extremely close wavelengths to Aq, namely A — 1499.9 
and 1500.2 nm. These correspond to L = 303.074 and 
303.133 /Ltm, respectively. 

Next, we consider the more realistic situation that L 
is fixed and the gain coefficient is adjustable. This is 
simply done by changing the pump intensity. It is not 



difficult to see that spectral singularities appear for spe- 
cific discrete values of go and A. This may be viewed as 
means for producing a tunable laser that would function 
at the very threshold gain. To implement this idea we 
insert ((23)) in ([20)1 and use ((22|) and A — ^-ncju) to express 
8r\ and k as functions of go and A. Substituting the re- 
sulting expressions in (fT7|) and (jTH]) and using fc — 27r/A 
we obtain a pair of equations for g^ and A for each value 
of the mode number m. The solution of these equations 
yield the desired values of go and A for which a spectral 
singularity appears. 

As we see from Figure [H the relevant range of values 
of br] and k is several orders of magnitude smaller than 
no. This shows that we can obtain reliable approximate 
forms of P7)l and (|18p by neglecting quadratic and higher 
order terms in br\ and k. Applying this approximation to 
([22)1 and using we have 



where 



/i := 



(1 -t2;2)2 +^^2^2 



h ■■= 



(1 -tj2)2 +^2^2 



(25) 



(26) 



w :— Lo/ujQ — Ao/A and 7 "j/i-Oo. In view of (|25]) . the 
above-mentioned approximation scheme is equivalent to 
neglecting terms of order two and higher in ko. Insert- 
ing ([25]) in ([T5)) and using k = 27r/A, we find 



Ko~ — 



h 



nLu}f2 
Ao 



(27) 



where po := ln[(?io + l)/(no - 1)]. In light of ^ 
is equivalent to 



go « 27rpo 



Ao/i 



1 



(28) 



Next, we substitute (1^ in ([T^ . neglect the quadratic 
and higher order terms in kq, and use (|27)) and p6)) to 
simplify the result. This yields after some lucky cance- 
lations: 



27rLnn 



m > u :~ 

' TT-V 



A, 

_ 2po 



Po 



TTmUJ 



Po 



0. Only for 
^ does this equation have real 



' 27rLytQ 
""7 V Po^o ^' 

solutions. These give the frequency (wavelength) of the 
spectral singularities. Substituting them in (j28l) . we find 
the corresponding go values. The approximate values of 
A and go obtained in this way allow for a more effective 
numerical solution of the exact equations (fT7|) and ([T8| . 

As a concrete example, we take L — 300 fiia for the 
gain medium described by (j24p and obtain the wave- 
length of the spectral singularities that are produced as 
we change go in the range 0-115 cm~^ . These turn out to 
correspond to the mode numbers 1355 - 1365 which are 
consistent with the result, m> p, k, 1320, of our approxi- 
mate calculations. The corresponding A and go values are 
given in Table H] The spectral singularity with lowest 
go value (40.4 cm~^) has a wavelength that is extremely 
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m 


A(nm) 


50 (cm ^) 


m 


A(nm) 


go (cm ^) 


1355 


1519.832 


110.2 


1361 


1495.662 


43.8 


1356 


1515.380 


82.5 


1362 


1494.582 


45.7 


1357 


1512.058 


66.3 


1363 


1489.203 


61.5 


1358 


1507.651 


50.9 


1364 


1484.927 


81.7 


1359 


1504.363 


43.8 


1365 


1481.736 


101.1 


1360 


1500.000 


40.4 









TABLE I: Wave length A of the spectral singularities for the 
gain coefficient go ranging over 0-115 cm~^ and L = 300 /im. 



sides of this equation gives rise to an independent con- 
dition for the existence of the spectral singularities that 
involves an integer (mode) number m. ft turns out that 
this equation and the threshold condition can be satisfied 
only for particular values of the physical parameters of 
the system and this corresponds to a single value of m 
and a corresponding critical wavelength. We have also 
explored the idea of tuning this wavelength by adjust- 
ing the pump intensity for a typical semiconductor gain 
medium. A remarkable feature of the spectral singularity 
related resonance effect is that if we increase the pump 



close to Aq = 1500 nm. It corresponds to m = 1360. As 
we increase go form 40.4 to 115 cm~^, there appear 10 
more spectral singularities with wavelengths ranging be- 
tween 1482 and 1520 nm. These should in principle be 
detectable, if we gradually increase the pump intensity. 
For example a periodic change of go in the range 70 - 
90 cm~^ should produce periodic emissions of radiation 
at the wavelengths A = 1484.927 nm (for m = 1364 and 
go = 81.6821 cm"!) and A = 1515.380 nm (for m = 1356 
and go = 82.4918 cm~^) with no emitted wave of com- 
parable amplitude at the resonance wavelength 1500 nm. 
This is a remarkable feature of the spectral singularity 
related resonance effect. Using the above values of A 
and go to compute the reflection and transmission coef- 
ficients and |Tp, that give the amplification factor 
for the emitted electromagnetic energy density, we obtain 
« |T|2 Tn 1.1 X 10^ and 3.1 x 10^ respectively. In 
other words we obtain an amplification of the background 
electromagnetic energy density at these wavelengths by a 
factor of |i?|2 + |T|2 « 2.2x10^ and6.2xl0^ respectively. 
It turns out that these numbers are extremely sensitive 
to the value of A but not go. Using the less accurate val- 
ues 81.7 and 82.5 cm~^ for go and the same values for A, 
we find sa |T|2 w 1.0 x 10*^ and 2.8 x 10^ But, as 
we can see from Figure [2l changing the above values of A 
by 0.01 nm reduces \R\^ and |Tp by three to four orders 
of magnitude. This shows that the detection of spectral 
singularities would require a spectrometer with a band 
width of 0.01 nm or smaller. 

To summarize, we have explored spectral singularities 
of a simple optical system and shown that the equation 
PU)) that determines spectral singularities reduces to the 
laser threshold condition g = gth, if we take the modulus 
of both sides of this equation. Equating the phase of both 
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FIG. 2: (Color online) Logarithmic plots of the reflection 
(thick dashed curves) and transmission (thin full curves) co- 
efficients as functions of A for L = 300 /im and g — 81.7/cm 
(top graph) and g = 82.5/cm (bottom graph). 



intensity so that the gain exceeds the threshold value, 
this effect disappears. This marks a clear distinction be- 
tween the zero-width resonances associated with spectral 
singularities and the usual resonances that we encounter 
in optical resonators. 
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